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A Burnett apparatus has been developed for the determination of the compressibility factor of 
gases, including corrosive halogenic gases, to an accuracy of 0.1 percent to 0.01 percent at temper- 
atures from to 225 °C and for pressures up to 250 bars. The apparatus has been used throughout this 
temperature and pressure range and at temperatures within 15 °C of the critical temperature of the 
sample gas. The apparatus features a rugged yet highly sensitive and reproducible null-type pressure 
transducer to isolate the sample gas without introducing a pressure uncertainty greater than the pre- 
cision of the pressure measurements, which is better than 1 part in 50,000. The pressure measurements 
are accurate to at least 1 part in 20,000. The reduction of the data is discussed extensively. The ca- 
pability of the apparatus has been evaluated in terms of the compressibility factor represented by a finite 
pressure or density virial expansion and derived from argon data at 25 °C and pressures up to 250 bars 
and in terms of the second virial coefficient of carbon dioxide based on data at 50 °C and pressures 
below 35 bars. The value obtained for this second virial coefficient, — 102.2 ±0.2 cm 3 /mol, compares 
favorably with the value of —102.5 cm 3 /mol obtained from a recent correlation and is in substantial 
disagreement with the value of —100.7 ±0.4 cm 3 /mol determined by Dadson and co-workers from the 
Burnett method and the piezometer method. 
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1. Introduction 

The accurate determination of the equation of state 
of a gas over a wide range of temperature and pressure 
provides important information about the gas both 
from a fundamental and an applied point of view. 
For example, the thermodynamic properties of a gas 
can be reliably extended beyond its ideal gas proper- 
ties if its equation of state is known with sufficient 
accuracy. Equation of state data are conveniently 
represented by means of the compressibility factor. 
The compressibility factor Z is defined in terms of the 
gas constant R and the measurable gas variables pres- 
sure P, volume V, absolute temperature 7 1 , and quan- 
tity, expressed as the number of moles n, by the rela- 
tionship 



Z = PV/nRT. 



(1) 



The accuracy of Z obtained from a straightforward 
measurement of each variable is limited by the 
accuracy with which the difficult measurements of 
n and V can be made. The direct or indirect measure- 



ment of one (n) or both (n and V) of the difficult vari- 
ables can be eliminated by the use of the virial expan- 
sions of Z as a function of P or of the density p and the 
fact that the dependence of Z on P is dominantly 
linear at sufficiently low pressures and approaches the 
ideal gas law 



Z = PV/nRT=l 



(2) 



in the limit as P goes to zero [1,2].* The Burnett PVT 
method [3-7] is an experimentally simplified method 
in which neither n nor Vis measured; instead, accurate 
isothermal pressure measurements are made before 
and after stepwise expansions between two volumes. 
The decrease in the overall measurement error 
achieved by the elimination of n and V measurements 
can in principle lead to more accurate values of Z 
than can be obtained from a method requiring the meas- 
urement of each variable in eq (1). For this reason we 
have developed a Burnett apparatus for the accurate 
determination of the compressibility factor of gases, 



1 Figures in brackets indicate the literature references at the end of this paper. 
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including corrosive halogenic gases [8], at tempera- 
tures from to 225 °C and for pressures up to 250 bars 
(1 bar= 10 5 N • m -2 ). The apparatus has been used by 
us throughout this temperature and pressure range 
including temperatures within 15 °C of the critical 
temperature of the sample gas. The apparatus is 
equally well-suited and has been used for the accurate 
determination of second virial coefficients. 

This paper is concerned with both the apparatus 
for obtaining Burnett data and the methods for re- 
ducing the data. These subjects are described briefly 
in the remainder of this section. In the next section the 
apparatus is discussed in more detail, and its dis- 
tinctive features are emphasized. Since the reduction 
of the data obtained with the apparatus is an essential 
element in the evaluation of its capability, the data 
reduction is discussed extensively in the following 
section. The discussion is illustrated with results for 
the compressibility factor of argon at 25 °C and pres- 
sures up to 250 bars and for the second virial coef- 
ficient of carbon dioxide at 50 °C. 

The apparatus, shown schematically in figure 1, 
consists essentially of a sample volume V x and an ex- 
pansion volume V u as well as components for (1) 
isolating, evacuating, or filling both volumes, (2) 
measuring the pressure in the sample volume, and (3) 
thermostating and measuring the temperature of the 
apparatus. The gas in the sample volume is isolated 
from the expansion volume by an expansion valve and 
from the external system by a metal diaphragm which 
is the pressure-responsive element of a null-type 
transducer. A counterbalancing gas acts on the other 
side of the diaphragm; the diaphragm is in its null 
position when the pressure on each side of the dia- 



phragm is the same. The pressure of the sample gas 
is obtained indirectly from the measurement of the 
pressure of the counterbalancing gas with a piston gage 
and a barometer. The sample and the expansion vol- 
umes include the appropriate volumes of the trans- 
ducer, expansion and filling-evacuation valves, and the 
interconnecting pressure lines; in our apparatus all 
these volumes are thermostated at a temperature which 
is regulated to about 0.001 °C and measured with a 
platinum resistance thermometer. Because no correc- 
tion needs to be made for temperature differences 
among these volumes, a correction that is often a 
significant source of error in PVT measurements is 
thereby eliminated. The transducer produces negli- 
gible variations of the sample volume and has a negli- 
gible effect on the high precision and accuracy of the 
pressure measurements. This is due to the excellent 
reproducibility of the null position of the diaphragm 
and the high sensitivity of the transducer. 

The experimental procedure is begun by filling 
the sample volume to the initial pressure, thermostat- 
ing the apparatus at the isotherm temperature T c , 
isolating the volumes from each other, and then 
measuring the gas pressure Po in the sample volume. 
The corresponding compressibility factor Zo is 



Zo = PoV l (Po)ln RT c , 



(3) 



where no represents the initial number of moles and 
Vi(Po) indicates that the sample volume is pressure 
dependent as a result of its distortion under pressure. 
In the next step, which we call the expansion step, 
the expansion volume is evacuated and isolated, 
and the gas in the sample volume is allowed to expand 
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FIGURE 1. Schematic diagram of the Burnett PVT apparatus. 

166 



into it. After temperature and pressure equilibration 
has occurred, the pressure P\ of the gas that now 
occupies both volumes is measured. For this first 
expansion step the compressibility factor Z\ is 



Zi=P x [V l (Px) + V ll (Px)VnoRTci 



(4) 



where /i and T c are the same as in eq (3). By taking the 
ratio of eqs (3) and (4) and defining N\ as the ratio of 
the combined sample and expansion volumes to the 
sample volume, i.e., N x = [F,(Pi) + fn(^i)]/^(^o), 
one obtains 



P n jP x = N,Z„IZ u 



(5) 



in which N\ depends slightly on both P () and P\. 
The ratio of the volumes can be written as 
Ni = N[l + bN^Po, Pi)], where N, known as the cell 
constant, is the value of the volume ratio in the limit as 
the pressure goes to zero and AA^i(P , P\) is a cell 
constant correction for the dependence of the cell 
volumes on pressure. 

Generally the sample volume is isolated again and 
the expansion step is repeated for the reduced quan- 
tity of gas rii that is now in the sample volume. This 
leads to a new pressure ratio, corresponding to eq (S), 



PilP2=N 2 Z 1 IZ 2 , 



(6) 



where N 2 = N[l+bN 2 (Pu Pi)]- This procedure is 
usually repeated to successively lower pressures with 
the successively smaller quantities of gas that remain 
in the sample volume after each expansion until the 
lowest pressures that can be measured both accurately 
and precisely are reached. The generalized form of eqs 
^5) and (6) for the rth expansion is 



Pr-ilPr = N r Zr-ilZr. 



(7) 



The combination of eq (7) and the relationship 
P = pZRT, obtained from eq (1), leads to a correspond- 
ing Burnett equation in terms of the density 



Pr = Pr-liN r . 



(8) 



Equations (7) and (8) are basic equations in the re- 
duction of Burnett data and may be transformed into 
different relationships in which the compressibility 
factor is represented by a finite virial expansion in either 
pressure or density. Any one of the relationships may 
be applied to Burnett data to form an overdetermined 
set of equations for evaluating the virial coefficients 
and generally the cell constant by means of nonlinear 
statistics [9]. Regardless of the relationship used, the 
cell constant is not determined from volume measure- 
ments. In the data reduction the pressure corrections 
ATVi, A7V 2 , . . . are not considered as unknown quan- 
tities but are estimated from theoretical and empirical 
distortion equations for pressure vessels or are deter- 
mined from a separate distortion experiment. 



2. Apparatus and Experimental Procedure 

2.1. Burnett Vessel and Gas Manifold System 

An overall view of the apparatus is shown in figure 2. 
In the Burnett vessel, as shown in figure 3, the sample 
chamber, 0.022 m i.d. by 0.076 m long, and the expan- 
sion chamber, 0.022 m i.d. by 0.051 m long, are bored 
from opposite ends of 0.064 m o.d. nickel rod to within 
0.051 m of each other. Nonrotating end closures are 
used for coupling to the transducer and to the filling- 
evacuation line. The top gland nut is split to accommo- 
date the transducer coupling. The coupling to the 
expansion valve is made through 0.0032 m o.d. by 
0.00080 m i.d. stainless steel capillary tubing silver- 
brazed into the sample chamber inlet. The coupling 
from the expansion chamber is made through similar 
capillary tubing silver-brazed into a ^-connection which 
in turn is coupled to the expansion valve and to the 
filling-evacuation valve. A well for the platinum resist- 
ance thermometer is bored radially into the Burnett 
vessel between the sample and the expansion chambers. 

Since halogenic gases react readily with most met- 
als, particularly if a trace of moisture is present, the 
components exposed to the sample gas are fabricated 
from metals such as nickel and stainless steel that are 
resistant to corrosion. We have found that gold provides 
still better resistance to corrosion than nickel and 
stainless steel do, and have therefore gold-plated every 
surface exposed to the sample gas in the Burnett vol- 
umes except the surfaces in the Burnett valves. Ex- 
treme care is taken to keep the system free of moisture 
at all times. The sample gas is not reclaimed; if neces- 
sary, neutralizing columns are used for its disposal. 
Cold traps prevent corrosive gases from entering the 
vacuum system, which includes low and high vacuum 
pumps, a cold cathode vacuum gage, and several ther- 
mocouple vacuum gages. Commercial diaphragm com- 
pressors boost the gas pressure from the supply cylinder 
pressure to about 600 bars for storage in high pressure 
receivers. 
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FIGURE 2. Burnett FVT apparatus. 
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FIGURE 3. Cross section of the Burnett volumes and valves. 



Carefully selected commercial nonrotating-stem 
valves are used; our experience shows that for a gas 
system the sealing surfaces of rotating-stem valves 
tend to gall and consequently to leak. The condition 
of the expansion and the filling-evacuation valves is 
particularly important; minimum deformation of the 
valve seat, long periods of trouble-free use, and 
certainty of valve closure are essential characteristics. 
A small valve with a block of 316 stainless steel, 
ports to accommodate 0.0032 m o.d. capillary tubing, 
and a stem of 17-4 stainless steel hardened to 40 
R/C is used for each of these two valves. Valve blocks 
for all the valves are selected on the basis of the 
concentricity of the seat and the stem-packing section. 
The shoulder of the valve seat is beveled slightly 
and carefully lapped to the same angle as that of the 
conical end of the valve stem. The sealing surface is 
kept small to minimize the torque required for valve 
closure. To prolong valve life, all valves are closed 
with a torque wrench, and the required torque is 
recorded on each valve block. 



2.2. Pressure Measurement System 

a. Transducer 

A rugged yet highly sensitive and reproducible 
null-type pressure transducer separates the sample 
gas from the pressure measurement system. The 
transducer is shown in figure 4; an earlier model has 
been described previously [10]. The transducer 
consists basically of a prestressed metal diaphragm 
as the pressure responsive element, a capacitance 
sensor for detecting deflections of the diaphragm as 
small as 0.025 ^m or less, and diaphragm backing 
surfaces incorporated as part of both the base and 
the capacitance sensor. 

The diaphragm is plane and circular, 25 fxm thick, 
and prestressed to reduce inelastic behavior. The 
backing surfaces, spherically concave with a chord of 
0.0238 m that determines the effective diameter of 
the diaphragm and a sagitta of 25 /urn, are necessary 
because the stepwise expansions in the Burnett 
method result in temporary pressure imbalances 
across the diaphragm many times greater than the 
elastic range of a diaphragm with the required 
sensitivity. 

The earlier model has been modified to improve its 
performance, especially above 125 °C. A metallized 
ceramic disk replaces the former metal-epoxy capac- 
itance sensor assembly. This has enabled us to operate 
the apparatus at temperatures up to 225 °C, and the 
disk could even be used at temperatures as high as 
500 or 600 °C. The diaphragm is now fabricated from 
rolled nickel, instead of Monel, sheet stock; its thermal 
expansion thus closely matches that of the nickel 
base, and the diaphragm tension does not depend 
significantly on the temperature. A polytetrafluoro- 
ethylene ring, prestressed by the clamping of the top 
and the base of the transducer and supported by 
inner and outer metal extrusion rings, replaces the 
former F-ring as the outer pressure seal. The dia- 
phragm sealing ring has been redesigned so that it 
stretches the diaphragm radially as it forms a pressure 
seal and thus does not distort the diaphragm when 
the transducer is assembled. Consequently the sample 
volume is more precisely defined and the transducer 
behaves more reproducibly. An important source of 
error is thereby reduced significantly. 

The transducer capacitance, consisting of the capac- 
itance contributed by the capacitor plate, diaphragm, 
and intervening gas plus the residual capacitance, 
constitutes one arm of a capacitance bridge. For 
each temperature and pressure the bridge is initially 
adjusted so that it is balanced when the pressure on 
each side of the diaphragm is the same. A subsequent 
pressure imbalance deflects the diaphragm from its 
null position and causes a bridge imbalance. The 
corresponding signal is amplified, phase detected, 
and indicated on a 100-0-100 /jlA output meter. 
The counterbalancing gas pressure is then adjusted 
until the diaphragm returns to its null position as 
indicated by the bridge balance. 
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FIGURE 4. Cross section of null -type pressure transducer. 

A. Follower, steel 

B. Extrusion rings, lava stone 

C. Top, nickel 

D. Clamping bolts for pressure seal K 

E. King, ceramic 

F. Capacitor lead stop, silver-brazed bead 

G. Clamping ring lor capacitance disk 

H. Inner clamping screws for clamping ring G 

I. Outer clamping screws for clamping ring J 

J. Clamping ring 

K. Polytetrafluoroethylene ring witb brass 

L. Capacitance disk 

M. Capacitor plate, metallized 

N. Diaphragm, nickel 

O. Capacitor lead, drill rod 

P. Electrical and oil sbield, Monel 

Q. Counterbalancing gas inlet tube, nickel 

R. Packing gland nut, steel 

S. Electrical insulating sleeve, polytetraflui 

T. Extrusion rings, ceramic 

U. Packing rings, polytetrafluoroethylene 

V. Metallized surface to capacitor plate 

W. Thrust screws for clamping capacitance disk 

X. Thrust screws for diaphragm pressure seal A- A 

Y. Thrust screws for stressing diaphragm N 

Z. Stressing ring, nickel 

A-A. Diaphragm seal 

B-B. Diaphragm stressing-sealing ring 

C-C. Diaphragm hacking surfaces 

D-D. Base 

E-E. Sample gas port 
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The temperature dependence of the transducer 
null capacitance is negligible for temperature changes 
even two orders of magnitude greater than the fluc- 
tuations of 0.001 °C about the isotherm temperature 
that represent the resolution of the temperature con- 
trol system. For each isotherm the transducer is 
calibrated before Burnett measurements are made: 
The counterbalancing gas, which has been argon in 
all of our work, is admitted to both sides of the dia- 
phragm, and the null capacitance is measured as a 
function of pressure from atmospheric pressure up to 
the highest Burnett pressure to be measured. The 
pressure dependence found does not represent any 
significant change in the diaphragm null position but 
is caused primarily by the change in the dielectric 
constant of the counterbalancing gas. The pressure 
sensitivity of the transducer, as used by us in the 
Burnett apparatus for pressures up to 250 bars, ex- 
ceeds the sensitivity of the piston gages and the 
barometer used to measure the pressure. 

b. Piston Gages and Barometer 

Although Bourdon gages are used, for convenience, 
as pressure indicators at various locations in the ap- 
paratus, the basic pressure measurements are made 
with piston gages and a standard barometer — instru- 
ments which are primary pressure standards. Two 
controlled-clearance piston gages with an accuracy 
of better than 1 part in 20,000 cover the range from 
about 2 to 3 bars gage pressure up to 250 bars and 
beyond. A gage with gas as the pressure fluid covers 
the lower part of the range, and a gage with oil as the 
pressure fluid covers the higher part. Both gages cover 
the region from about 30 to 45 bars and in this overlap 
region they agree with each other to within 1 part in 
40,000 to 50,000 of the measured pressure. 

The nominal effective area of the oil gage piston is 
26 X 10 -6 m 2 and of the gas gage piston, about 
340 X10 -6 m 2 . The pistons are fabricated from tungsten 
carbide; the cylinders, from stainless steel. The top 
of the oil gage piston is hemispherical and fits into a 
conical recess in the coupling that connects the piston 
to the weight holder on which the weights are loaded. 
In the gas gage a lubricated ball bearing fits into a 
conical recess in both the piston and the coupling, in 
each case the purpose of the spherical contact, which 
represents our modification of the gages, is to reduce 
the amount of side thrust applied to the piston. Simi- 
larly, the weight holders and the weights that nest on 
them are machined carefully to avoid eccentric loading. 
The prevention of side thrust is particularly important 
for the gas gage because its piston receives little lateral 
support. The gas gage piston and cylinder must also 
be kept scrupulously clean. For each gage the load 
applied to the piston is spun manually and continues 
to spin for up to an hour or more before it stops. 

Even though the pistons and cylinders were care- 
fully lapped and cleaned before they were delivered 
to us, we were able to improve the performance of the 
gages by using a simple motor drive assembly to rotate 
the piston and move it up and down within its cylinder 
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for several hours while the piston and cylinder were 
flooded with mineral spirits. The resultant burnishing 
increased the spin time, and the flooding helped to 
remove any fine particles that might have been present. 
The burnishing procedure can be used to improve 
performance at any time, particularly if the cylinder 
has become slightly scratched. 

A simple piston gage without the controlled-clear- 
ance capability is used to make less accurate meas- 
urements over the range from 2 to 3 bars gage pressure 
down to 0.2 bars absolute pressure or less; the abso- 
lute measurements are made under an evacuated 
bell jar. The gage has three stainless steel piston and 
cylinder assemblies with nominal effective areas of 
800 X 10- 6 m 2 , 80 X 10" 6 m 2 , and 8 X 10" 6 m 2 , and uses 
gas as the pressure fluid. 

The atmospheric pressure is measured with a stand- 
ard barometer of modified U-tube design consisting of 
a column and a cistern. The level of the mercury in the 
column is detected photoelectrically and is measured 
with a 31 in (0.79 m) scale and a vernier graduated to 
0.001 in (25 /xm); the level of the mercury in the cistern 
is detected electrically and is measured with a microm- 
eter graduated to 0.0001 in (2.5 fjon). The back pressure 
above the mercury in the column is measured with a 
McLeod gage and is generally about 4 X 10 -6 bars or 
less. The accuracy of the barometer scale is easily 
checked with gage blocks. Repeated tests have shown 
that the barometric pressure measurements are ac- 
curate to 33 X 10~ 6 bars. 

c. Other Components 

Small adjustments in the counterbalancing gas pres- 
sure are made with a compound gas injector whose 



internal volume is changed by the manual displace- 
ment of a 0.013 m or 0.0032 m diam piston. The coun- 
terbalancing pressure acts directly on the gas gages but 
is separated from the oil gage either by a commercial 
diaphragm transducer or by a mercury U-tube. The 
U-tube is an adaptation of an earlier model by Keyes 
[11] and is shown in figure 5. To reduce capillary effects, 
the lower part of each bore containing the mercury is 
made large in diameter. The bore diameter is reduced 
slightly above the mercury to make the gas volume in 
the bore a small fraction of the total volume of the 
counterbalancing gas. As a result, the bore gas volume 
has only a small effect on the performance of the gas 
injector. The mercury level in each column is indicated 
by electrical contact between the mercury and a fine 
metal probe whose height is initially adjusted to within 
25 /mi of the reference height and is checked period- 
ically. The modified design permits the valve at the 
base of the U-tube to be replaced conveniently. 

2.3. Temperature Control and Measurement System 

An oil bath and a servo heater system are used to 
control the temperature of the Burnett apparatus. As 
shown in figure 6, the Burnett vessel, the transducer, 
the expansion and the filling-evacuation valves, and the 
interconnecting lines are all immersed in oil in an inner 
chamber located in a thermally insulated bath. The 
inner chamber is double-walled, and the space between 
the walls can be evacuated to reduce radial heat flow. 
An oil-circulating pump at the bottom of the bath is 
coupled mechanically to a motor mounted on the top. 
Oil drawn into the manifold and the pump is ejected 
upward into the inner chamber and flows by the Burnett 
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Figure 5. Mercury U- tube. 
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Figure 6. Oil bath. 

apparatus, over the top of the inner chamber, and then 
downward past the heat exchanger to complete the 
cycle. Sufficient mixing of the oil occurs in the mani- 
fold and the pump to eliminate temperature gradients 
acquired during the flow of the oil past the heat 
exchanger. 

The heat exchanger consists of a cooling coil and a 
series-parallel heater for low and high power output. 
A coolant is circulated through the coil and a refrig- 
erated bath, which is thermostated to within 0.003 °C 
from —20 to 5 °C. The cooling system is used for oil 
bath temperatures from °C to slightly above room 
temperature. 

The servo heater system consists of an adjustable 
Wheatstone bridge with a resistance thermometer con- 
stituting one arm, a d-c amplifier and recorder, a three- 
mode controller for proportional, derivative, and inte- 
gral response to the error signal, and a power amplifier 
connected to a heater circuit. The resolution of the tem- 
perature control is 0.001 °C; slight differences from the 
desired temperature are readily corrected by the use 
of a variable resistor in the bridge circuit. 

The temperature as defined by the International 
Practical Temperature Scale of 1968 is measured with 
a capsule platinum resistance thermometer located in 
the well between the sample and the expansion cham- 
bers in the Burnett vessel. The thermometer leads are 
brought out through a protection tube filled with oil to 
ensure good thermal contact between the thermometer 
and the Burnett vessel. A Mueller bridge with a photo- 
electric galvanometer as the detector is used to measure 
the resistance. 



Differential five-junction copper-constantan thermo- 
couples are used to measure temperature differences 
over the Burnett vessel and transducer as well as near 
the top of the oil bath. A small temperature difference 
can exist between the Burnett vessel and the top of 
the bath due to heat loss through the apparatus supports 
and oil-vapor loss to the surroundings. This difference 
is limited to 0.01 °C by means of a manually-controlled 
auxiliary heater at the top of the bath. Temperature 
differences over the Burnett vessel are less than 
0.002 °C. Electrical noise in the thermocouple circuits, 
which include the switching and measurement instru- 
mentation, corresponds to a temperature uncertainty 
of less than 0.0001 °C. 

2.4. Experimental Procedure 

An experimental procedure has been evolved which 
ensures the reproducibility of the pressure measure- 
ments, the precise definition of the sample volume, 
and the detection of any sample gas leakage. A series 
of runs is made at each isotherm temperature. Each 
run may include up to approximately eight expansions; 
the initial pressure for each run is chosen so as to 
give a well-distributed set of pressures for the iso- 
therm. Before any Burnett measurements are made, 
the oil bath is thermostated at the desired tempera- 
ture and the transducer is calibrated over the desired 
pressure range. The transducer null capacitance at 
atmospheric pressure is checked before and after 
each day's measurements. 

The first step in a run is to evacuate the sample gas 
part of the apparatus and flush it thoroughly with the 
sample gas. After this part of the apparatus has been 
reevacuated, the sample and the expansion volumes 
are filled with fresh sample gas to the highest pressure 
desired for the run. In the flushing and the filling proc- 
ess and throughout the run, any pressure imbalance 
across the transducer diaphragm greater than a few 
bars is directed toward the base, which is more rugged 
than the capacitance sensor and its clamping arrange- 
ment and hence provides better diaphragm support. 
This precaution helps to ensure the stability of the 
diaphragm null position and thus aids in the precise 
definition of the sample volume. Similarly, as a pres- 
sure measurement is made, the diaphragm is slightly 
overloaded toward the base prior to the final balancing 
of the counterbalancing and the sample gas pressures. 

After the filling process is completed and the filling- 
evacuation valve closed, sufficient time is allowed 
for temperature and pressure equilibration. During 
the course of a run the temperature of the Burnett 
vessel may vary over a range of about 0.002 °C due to 
instabilities primarily in the temperature control 
bridge; equilibration times are chosen to keep short- 
term temperature changes well within this range and 
temperature differences within I he ranges discussed 
earlier. 

While the equilibration is taking place, preparations 
are made for the first set of Burnett measurements. 
The counterbalancing gas, which was maintained at 
a higher pressure than the sample gas during the filling 
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process, is gradually vented until the transducer null 
indicator returns to its zero position. The approximate 
counterbalancing gas pressure is then read from a 
Bourdon gage, and a variable capacitor in the capaci- 
tance bridge is adjusted in accordance with the 
transducer calibration so that the bridge balance will 
truly correspond to the diaphragm null position. The 
piston gage to be used for the precise pressure meas- 
urement is loaded with approximately the right 
combination of weights that will be needed. 

When the resistance thermometer and the differen- 
tial thermocouple readings show that the equilibration 
is essentially complete, the expansion valve is closed 
to isolate the sample and the expansion volumes. While 
the valve is slowly closed, the diaphragm is maintained 
in its null position by adjustment of the counterbalanc- 
ing gas pressure with the compound gas injector. An 
amplifier and supplementary null indicator about ten 
times as sensitive as the transducer output meter are 
used in this process, partly because of their greater 
sensitivity and also because the location of the supple- 
mentary null indicator is more convenient to the person 
closing the valve and operating the compound gas 
injector. 

After the expansion valve is closed, the first set of 
Burnett measurements can be made. In the operation 
of the piston gage a precise zeroing of the transducer 
output meter would require an unnecessarily long time 
to be achieved. Instead the gage operator calibrates 
the output meter over its central range by noting the 
difference in the scale readings that correspond to two 
nearly equal loads on the gage. The slight difference in 
the diaphragm position represented by the different 
scale readings changes the sample volume, which is 
bounded by the diaphragm, by a negligible amount — 
less than 2 ppm. The scale readings are made when 
the piston is at its reference height as shown by an elec- 
tronic proximity sensor mounted underneath the weights 
loaded on the gage. The height of a gas gage piston is 
adjusted by means of the compound gas injector; the 
height of the oil gage piston is adjusted by means of a 
separate oil injector. For either type of gage the piston 
is raised above its reference height by means of the ap- 
propriate injector and is then allowed to fall slowly down 
to its reference height while the piston and its load are 
spinning steadily. The gage pressure and barometric 
pressure measurements are made almost simultane- 
ously, and the resistance thermometer and differential 
thermocouple measurements are made immediately 
afterward. These measurements constitute the first 
set of Burnett measurements. 

Next the expansion valve is opened and after about 
a minute is closed as before. A second set of Burnett 
measurements is made. Since the expansion step has 
not been performed, this second set of measurements 
should be in close agreement with the first set. If the 
second set, in fact, does not differ significantly from the 
first, the filling-evacuation valve is opened and the ex- 
pansion volume is evacuated. After isothermal condi- 
tions are reestablished and a thermocouple vacuum 
gage, located as close to the expansion volume as is 
convenient, shows a pressure of about 4 X10~ 6 bars, a 



third set of measurements is made. If this set does not 
differ significantly from the first two, the filling- 
evacuation valve is closed and the expansion valve is 
opened. The gas initially in the sample volume ex- 
pands to fill the expansion volume as well. After this 
first expansion the measurement procedure is repeated 
to successively lower pressures. 

To be considered satisfactory, the three sets of meas- 
urements corresponding to the initial filling or a sub- 
sequent expansion must show a pressure reproducibility 
of 1 part in 50,000 or better. Occasionally a set of meas- 
urements will be repeated if some temporary disturb- 
ance, such as a temperature or barometric fluctuation, 
has caused a poorer reproducibility. If a serious ex- 
perimental failure, such as sample gas leakage, were 
to occur, the measurements for that run would not be 
used in the data reduction. 

3. Data Reduction and Discussion of 
Results 

3.1. Background 

In contrast to the simplicity of the experimental 
measurements, the reduction of Burnett data involves 
rather formidable nonlinear statistical analysis of 
relationships containing the compressibility factor, 
defined by a pressure or a density virial expansion, 
and optional constants. The authors have previously 
presented a general discussion of the analysis of 
different Burnett relationships using two nonlinear 
statistical methods, due to Gauss and to Deming [9], 
Related papers have also been published by others 
[12, 13]. 

Our experience indicates that compressibility 
factors and second virial coefficients calculated from 
Burnett data of high precision and accuracy are as 
accurate as those calculated from the best data ob- 
tained with any other PVT method. Various treatments 
of relationships, that is, various models, based on eqs 
(7) and (8) have been used in our analysis. In these 
models the cell constant, where it appears, is regarded 
as a single parameter common to all the experimental 
runs, as a separate parameter for each run, or as a 
constant to be determined from Burnett measurements 
for a nearly ideal gas. Meaningful comparisons of 
the results are difficult for at least two reasons: (1) 
the same systematic errors may propagate differently in 
each model and (2) the number of parameters to be 
evaluated varies from one model to another. The 
effect of the latter is subtle; it may introduce an 
analytical bias, as the magnitude of the correlation 
between the parameters, which is basically high, 
increases significantly with an increase in the number 
of parameters. It is not surprising, therefore, that slight 
discrepancies exist among the results of different 
models, even when applied to the same data. For 
example, a comparison of the compressibility factors 
derived from argon data taken with our Burnett 
apparatus shows discrepancies among the values of 
approximately 0.01 percent to 0.04 percent, which, 
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although small, are nevertheless greater than the 
estimated statistieal uneertainties of less than 0.01 
percent. For the six models considered, the results of 
four show comparative discrepancies of about 0.01 
percent. Another example of this type of discrepancy 
occurred in an informal round-robin analysis of helium 
data by participants at a two-day Burnett conference 
held at the National Bureau of Standards in February 
1969. Initial analyses indicated discrepancies of 0.04 
percent among the results of the different models as 
compared to statistical uncertainties of 0.01 percent 
in the results themselves. 

We shall not repeat our previous detailed presenta- 
tion of the nonlinear statistical analysis of Burnett 
data [9], but shall present a brief discussion which 
emphasizes the important differences between the 
various models in order to help one choose the models 
that introduce the least bias. The fact that the data can 
be reduced in several ways is an attractive feature of 
the Burnett method and enables one to study the dif- 
ferent effects of systematic errors. 

The basic Burnett equations may be transformed into 
different relationships among the variables of interest. 
In terms of a finite virial expansion in pressure, 

Z= 1 + V biP*, eq (7) may be expressed as 
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In terms of a finite virial expansion in density, 

j 
Z = 1 + V B,pJ, eq (8) may be expressed as 
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It should be noted that in these equations the second 
virial coefficient is represented by b\ or B\. For both 



the pressure and the density relationships 
Nr = N[l+bN r (Pr-uPr)] 

and 

NiN, . . . /V r =/V'[l+A,V, (Po, Pi)] 

X [l + AiV|(Pi,Pt)] ••• [I +AN r (Pr-l, Pr)] . 

Any one of these relationships may be applied to 
Burnett data to form an overdetermined set of equa- 
tions for evaluating the virial coefficients and the 
other parameters. The latter consist of the cell con- 
stant (except in eq (11) where the cell constant has 
been eliminated), the initial density po in eq (12), 
and sometimes the quantity Po/Zq in eq (10). The latter 
two parameters are characteristic of an individual 
run. As mentioned in the introduction, the pressure 
corrections (ATVi, A/V 2 , . . .) are not considered as 
parameters but are determined independently of the 
parameter evaluation. 

3.2. Compressibility Factor of Argon at 25 C 

Compressibility factor results derived from the 
different relationships applied to our argon data are 
given in table 1 [9J. The models are designated briefly 
at the top of each column. For each model the param- 
eters were evaluated simultaneously for all the data. 
Where a fixed cell constant was used, it was pre- 
determined from eq (9) applied to helium data taken 
at the same temperature for pressures below 35 bars. 
For these conditions helium behaves nearly as an 
ideal gas, a characteristic which simplifies the data 
reduction. 

The results at each pressure are essentially the same 
except for those derived from the density relationship 
with an adjustable cell constant included as either 
a separate parameter for each run or a single param- 
eter common to all the runs. The discrepancies among 
the first four columns, which are about 0.01 percent, 
could be due to the presence ot systematic error which 
propagates differently in each model. A likely source 
of such error is the correction for the volume dis- 
tortion of the sample and expansion chambers with 
pressure. A study of simulated Burnett problems 
based on our argon data indicated that other probable 
sources of error had a negligible effect on the com- 
pressibility factor. In this study a simulated syste- 
matic error of 1 part in 30,000 was superimposed on 
the argon pressure data. The density reduction of 
these data, with the cell constant either fixed or re- 
garded as a parameter, resulted in a change in the 
compressibility factor of less than 0.001 percent for 
all the pressures. An effect of less than 0.001 percent 
was also obtained when a simulated systematic error 
of 1 part in 20,000 was superimposed on the RT product 
in a density reduction of the data. 

The volume distortion correction was determined 
from distortion equations and strain gage measure- 
ments of the external deformation of the Burnett 
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Table 1. Comparison of results derived from Burnett argon data at 25 °C and to 250 bars a ' h 



3 Runs (7, 7. & 6 expansions/run) 





Pressure series formulation 


Density series formulation 




Eq(9) 


Eq(9) 


Eq(ll) 


Eq (12) 


Eq (12) 


Eq (12) 




Cell constant 


One cell constant 


Cell constant 


Cell constant 


One cell constant 


All cell constants 




fixed 


adjusted 


eliminated 


fixed 


adjusted 


adjusted 


N 


1.74505 


1.74502 




1.74505 


1.74511 


1.74511 
1.74509 
1.74511 


^10.133 bars 


0.99384 


0.99387 


0.99388 


0.99380 


0.99375 


0.99376 


Zioi.33 bars 


.95356 


.95365 


.95367 


.95352 


.95333 


.95334 


^202.85 bars 


.95202 


.95209 


.95212 


.95194 


.95169 


.95173 



a The compressibility factor is shown to five digits to indicate the differences in the numerical results and not to indicate the overall 
accuracy. 

b 1 bar= 10 5 N-m" 2 . 



vessel due to internal pressure [14, 15]. The effect 
of the distortion correction and its estimated uncer- 
tainty on the compressibility factor is shown in figure 
7 for argon data analyzed with the use of the pressure 
equation (10). The effect of the uncertainty is small, 
about 0.01 percent at the maximum pressure of 250 
bars. The results of simulated problems indicate that 
the effect is the same in the density formulation. In 
the pressure series formulation the effect of the dis- 
tortion correction uncertainty is consistent with the 
small differences between the results derived with the 
cell constant fixed and those derived with the cell 
constant evaluated or eliminated. 

We now consider the basis for the differences, 
particularly at the higher pressures, which are sig- 
nificantly greater than 0.01 percent. These are between 
the results derived from the pressure series formula- 
tion and the results from the density series formula- 
tion, particularly with the cell constant included as 
either a separate parameter for each run or a single 
parameter common to all the runs. If the exact model 
were known, that is, if all the systematic errors 



5 — 



4 — 



< 
o 



3 — 



1 1 ' 1 


1 | 1 | i | i | . | 


— 


^''CORRECTION 


- 


/ ESTIMATED UNCERTAINTY 


— / 


IN THE 


- / 


CORRECTION _ 


^^-^^ 


> 1 , 1 l 1 , 1 1 1 



40 



80 



120 



160 
bars 



200 



240 



280 



FIGURE 7. Relative compressibility factor correction due to pres- 
sure distortion of Burnett cell volumes, argon run. 



were included in each equation, then the uncertainties 
in the values of the compressibility factor, the cell 
constant, and the initial densities would be consistent 
with the magnitude of any random errors in the pres- 
sure. This has been illustrated with a simulated Burnett 
problem for a particular model in which random errors 
were superimposed on the pressure data [16]. In 
the real world one cannot account exactly for all of 
the systematic errors. Consequently, a bias is intro- 
duced in the statistical analysis, and the results are 
not as reliable as the statistical uncertainties indicate. 
In addition, the number of parameters varies from one 
model to another. This variation adds to the bias since 
the absolute value of the correlation between the 
parameters increases significantly with an increase 
in the number of parameters. 

In general the correlation between virial coefficients 
is high regardless of which PVT method is used to 
determine them. In the Burnett method additional 
highly correlated parameters are also present. In the 
density reduction of the argon data with a fixed cell 
constant, the absolute value of the computed correla- 
tion between each of the virial coefficients is about 
0.85 on a normalized basis of to 1, whereas it is about 
0.6 between each of the initial densities and about 
0.75 between any one of the virial coefficients and any 
one of the initial densities. In contrast, if the cell 
constant is included as a parameter, the statistical 
fit is improved over that for the fixed cell constant, 
in the sense that the sum of squares is reduced. 
However, the correlation between any two of the 
parameters increases significantly! The correlation 
between each of the virial coefficients becomes about 
0.93; between each of the initial densities, about 
0.994; and between any one of the virial coefficients 
and any one of the initial densities, about 0.87. The 
correlation between the cell constant and any one of 
the virial coefficients is about 0.88 and between the 
cell constant and any one of the initial densities, about 
0.995. Exceptionally high correlations such as 0.994 
and 0.995 indicate that these parameters are not truly 
independent and thus cannot be accurately evaluated 
[17]. 
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Additional data, not highly correlated to the data 
already in the analysis, would reduce the high param- 
eter correlation. As a result, the evaluation of the 
parameters would be enhanced and any analytical 
bias in the results would he reduced. In one approach 
the analysis would consist of data from the same ap- 
paratus for two or more sample gases with different 
compressibility characteristics. The data would be 
taken at the same temperature, and a single cell 
constant parameter, common to all the data, would be 
evaluated. Another approach would be to use a fixed 
cell constant obtained from helium data. This approach 
is less desirable because the cell constant would be 
entered without any statistical uncertainty, that is, as an 
exact quantity, and hence the computed statistical 
uncertainties would be too small. The significance of 
high parameter correlation should be considered care- 
fully in estimating the reliability of results obtained 
from the density reduction methods. For the pressure 
reduction methods the situation is similar, although the 
parameter correlation is not as high. 

Our values of the compressibility factor shown in 
table 2 represent the average of the results of the three 
pressure models and the density model with a fixed 
cell constant; the results of the density models with one 
or three separate cell constants were not included in 
the average because they were considered to be less 
reliable due to the high correlation between the 
parameters. For pressures up to 250 bars the overall 
accuracy of the averaged compressibility factor is 
estimated to be within ±0.0002 on the basis of our error 
analysis. The statistical uncertainty of each value of the 
compressibility factor is no greater than one sixth 
of the overall accuracy of the averaged compressibility 
factor. Table 2 also includes comparative values of 
the argon compressibility factor derived from the 
piezometer results of Michels [18] for pressures up 
to 202.7 bars. The maximum difference is small, only 
0.05 percent at 202.7 bars; PVT results are generally 
considered to be accurate if the overall error is less 
than 0.1 percent. 

Table 2. Comparative values of the compressibility 
factor for argon at 25°C ii 



Pressure 
bars 


This paper h 
Compressibility factor 


Michels 1 " 
Compressibility factor 


10.133 
101.33 
202.65 


0.9938 
.9536 
.9520 


0.9938 
.9534 
.9515 



a 1 bar= 10* N- m 2 . 

b The overall accuracy is estimated to be ±0.0002. 
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3.3. Second Virial coefficient of Carbon Dioxide 
at 50 °C 

In addition to accurate values of the compressibil- 
ity factor, reliable second virial coefficients are also 
obtainable from the Burnett method. Our own ex- 
perience with argon at 25 °C and with carbon dioxide 



at 50 °C, 19 degrees above its critical temperature, 
shows that this is true over a wide range of reduced 
temperatures greater than unity. 2 The argon results 
are for a reduced temperature of 1.98 and have been 
discussed previously [9]; the carbon dioxide results, 
which are discussed in the following paragraphs, 
are new and are for a reduced temperature of 1.06. 

Two groups of carbon dioxide measurements were 
made: one group over the pressure range from 35 to 
4.5 bars and another group over the range from 4.5 
to 0.5 bars. The higher pressure data are the more 
accurate data and were reduced by the density model 
with a fixed cell constant determined from helium data 
at the same temperature. From the more accurate 
data a second virial coefficient for carbon dioxide 
of —102.25 ±0.04 cm'Vmol was obtained, where the 
uncertainty is only the statistical uncertainty and is 
necessarily too small because it does not include the 
uncertainty in the cell constant. Our estimate of the 
overall systematic uncertainty, based on our error 
analysis, is 0.2 cm'Vmol. This value of the second 
virial coefficient was supported by a graphical evalu- 
ation which included the less accurate carbon dioxide 
data for the lower pressures. It compares favorably 
with the value of —102.5 cm'Vmol obtained from a 
recent correlation [19] which included an evaluation 
of all existing data but is in substantial disagreement 
with the value of —100.7 ±0.4 cm 3 /mol determined 
by Dadson and his co-workers [20, 21]. It should be not- 
ed that Dadson et al. carried out their measurements 
with both the Burnett method and the piezometer 
method. The fact that their results from both these 
methods agreed with each other and yet disagreed 
with the results of the correlation prompted our own 
measurements. Our results and more detailed compari- 
sons will be the subject of a later publication. 

Various potential sources of systematic error were 
investigated in our own work. In simulated problems 
based on the higher pressure data, superimposed 
systematic errors of ±25 percent in the distortion 
correction and 1 part in 20,000 in the RT product had 
only a negligible effect on the second virial coefficient. 
The uncertainty in the distortion correction would be 
expected to have a negligible effect since the distortion 
correction itself is very small even at the highest 
pressures in the carbon dioxide data. A cell constant 
was computed for the carbon dioxide data in a separate 
calculation by treating it as an adjustable parameter 
in the density series formulation. The value obtained 
agreed with the predetermined value obtained from 
helium data to within 1 part in 20,000, which is excellent 
agreement in view of the large magnitude of the non- 
ideality of carbon dioxide at 50 °C (777\. ritk . al =1.06) and 
of the extremely high correlation between the param- 
eters in the analysis. 

The possibility of adsorption of carbon dioxide on the 
surfaces of the sample and the expansion volumes was 
investigated by means of a graph of the low pressure 
data in which the pressure ratio P r -ilP, was plotted 



2 The reduced temperature is the ratio of the sample gas temperature to its critical 
temperature, where both these temperatures are expressed on the absolute scale. 
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as a function of the pressure P r -i- The significance 
of the graph is based on the equation 



Pr- 



Pr 



••N+bt(N-\)Pr- U 



(13) 



which is the low pressure form of eq (9). 3 The extrapo- 
lation of the curve drawn through the points to the verti- 
cal axis yielded an intercept whose value, the cell con- 
stant, was essentially the same as the predetermined 
value obtained from the helium data. This agreement 
would not be expected if adsorption were present. 
For a nonpolar gas such as carbon dioxide and for our 
range of pressures, adsorption is expected to be sat- 
urated [22]; that is, above a relatively low pressure 
very little more adsorption occurs as the pressure is 
increased further. In the sample volume the adsorbed 
gas would remain essentially unchanged upon an ex- 
pansion, whereas in the expansion volume gas would be 
desorbed in the evacuation process and adsorbed 
in the expansion process. As a result, the number 
of free molecules would be reduced upon an ex- 
pansion, and the ratio P r -i/P r would therefore in- 
crease rapidly at low pressures. The extrapolation 
to zero pressure would yield a value for the cell 
constant that is greater than the value predetermined 
from helium data for which adsorption is not signifi- 
cant. Furthermore, the graph of this ratio versus pres- 
sure would show increasing nonlinearity at low pres- 
sures instead of the decreasing nonlinearity character- 
istic of real gases as the effect of the third and higher 
order virial coefficients not shown in eq (13) becomes 
insignificant. Neither of these signs of adsorption is 
present in our data. Hence we feel confident that any 
adsorption in our apparatus is so slight as to have an 
insignificant effect on our results. 



We gratefully acknowledge the benefit we have 
received from general discussions with M. Klein and 
from comments regarding adsorption made by P. T. 
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3 Equation (13), together with the equation B t = RTb^, was used in the graphical evaluation 
of the second virial coefficient. 
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